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Abstract: 

Solid state quantum condensates can differ from other condensates, such as Helium, ultra- 
cold atomic gases, and superconductors, in that the condensing quasiparticles have relatively 
short lifetimes, and so, as for lasers, external pumping is required to maintain a steady state. In 
this chapter we present a non-equilibrium path integral approach to condensation in a dissipative 
environment and apply it to microcavity polaritons, driven out of equilibrium by coupling to 
multiple baths, describing pumping and decay. Using this, we discuss the relation between 
non-equilibrium polariton condensation, lasing, and equilibrium condensation. 
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1 Introduction 

The idea of Bose-Einstein condensation of quasiparticles in solid-state structures has a long 
history, since the early proposals [IHS] that excitons might form a condensate. In more recent 
years, this has led to work on a variety of systems: excitons in coupled quantum wells [1H7]; 
excitons in quantum-Hall bilayers jHj; magnons, both in thermal equilibrium [9] and parametri- 
cally pumped magnetic insulators |10| as well as within superfluid 3 He|lH [T2]; and microcavity 
exciton-polaritons. (For extensive references to experiments see the review [H])- In almost all 
these cases, the condensate is, to a greater or lesser extent, a non-equilibrium steady state, with 
pumping compensating for the finite lifetime of the quasiparticles, leading to a flux of particles 
through the system. Thus, a general question arises: can Bose-Einstein condensation be realised 
in a strongly dissipative environment, and if so how would it relate to and differ from equilibrium 
BEC and the laser? 

To address these questions, this chapter discusses a field theoretical approach to modelling 
quantum condensates that are driven out of equilibrium by a flux of particles through the system. 
We illustrate the technique using an example of microcavity polaritons. Their part-light nature 
leads to a rather short lifetime that may nonetheless be long enough to have polaritons as 
well-defined quasiparticles (i.e. strong coupling). Their short lifetimes however lead to an 
important role of non-equilibrium physics. This naturally provokes questions about the relation 
to lasing, which occurs in pumped cavities in the weak coupling limit. We consider a polariton 
system coupled to baths which model pumping and decay processes. Since these baths are 
not in chemical equilibrium with each other, they drive a flux of particles. The Hamiltonian 
we use will describe both a laser (if pumped at high temperatures, as discussed below), and 
Bose condensation if treated in thermal equilibrium, as well as the smooth transition between 
them. As such, the system of microcavity polaritons provides a particularly rich playground for 
studying coherence in a dissipative environment, and exploring the differences and similarities 
between condensates and lasers. 

2 Methodology: Modelling the Non-Equilibrium System 

For a non-equilibrium system, the density of states and its occupation must both be determined 
explicitly, as the occupation may be non-thermal. This means that to describe the system 
fully, one needs at least two Green's functions. We choose here to work with the retarded 
and Keldysh Green's functions: D R (r, r', t, t') = -i0(i)([^(r, i), ^(r', i')]-}, D K (r, r', t, t') = 
— i(fy(r, t), ^(v', w h ere [Vs V^]=f ^ s the commutator (anti-commutator). The retarded 

Green's function describes the response following some applied perturbation. In the frequency 
domain, p(p,c<j) = 2lm{D R (p,Lu)] gives the density of states, while the Keldysh Green's function 
D K (p,u) = —i[2n(uj) + l]p(p,cj) accounts for occupation n{ui). 

To determine these Green's functions, we will use a path integral approach[14J, discussed 
further below. Path integrals naturally allow computation of time-ordered correlation functions; 
in order to instead find the retarded and Keldysh Green's functions we must use the Keldysh 
contour C. Points on this contour are labelled by (t, {+,—}), where +, — distinguishes the 
forward(+) and backward(— ) branches. The path integral approach will then give contour or- 
dered correlations, denoted by Tq, where fields on the + contour always precede those on the — 
contour, and fields on the — appear in time reversed order. Then introducing symmetric (classi- 
cal) and anti-symmetric (quantum) combinations of these fields V>ci,q = [VK+iO ^ V'( — >*)] /v2, 
the Green's functions are given by: 

D =(o* T)-K^(«M))(* r ' 0) '* r ' 0, )>- (1) 
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(D A is the advanced Green's function, the Hermitian conjugate of D R ). As discussed below, the 
action in the path integral involves the inverse Green's function: D -1 = {j^A D q^ = 

( [2 ,K]-i [£-ijic)> where t^" 1 ]^ = -[D R ]~ 1 D K [D A ]- 1 . As an illustration, for a free field 
[D^] -1 = hu} — hu! p + i5 and [-D^" 1 ]^ = 2i5[2nB(fioo) + l], where 5 is infinitesimal. The above is for 
bosonic fields; the results for fermionic fields are similar, but commutators and anti-commutators 
are interchanged in the definitions of Keldysh and retarded Green's functions. 



Polariton System Hamiltonian, and Coupling to Baths 

To describe the polariton system we use a model of disorder localised excitons strongly coupled 
to cavity photons |15H17j . Exciton-exciton interactions are included in this model by allowing 
only zero or one excitons on a given site, thus describing hard-core bosons. This model has 
several advantages for our aims: Firstly, this same Hamiltonian has been used to model lasers 
|18j . allowing us to relate polariton condensation to lasing. Secondly, it is known [16] that, in 
equilibrium, except at extremely low densities, mean-field theory captures the phase diagram of 
this model rather well. Finally, it allows one to account straightforwardly for exciton nonlinearity 
within the non-equilibrium mean-field theory. 

To describe a hard-core boson, we introduce two fermionic operators dt, ct that create states 
representing the presence or absence of an exciton. The operator djc- thus creates an exci- 
ton. In this notation, the system Hamiltonian is H sys = J2j e j($jdj ~ + Yip hup^p^p + 
Yj P 9j(^pCjdj + H.c.). Here ej is the exciton state energy and gj is the coupling to photons. 
The cavity photon dispersion is hco p = Hloq + h 2 p 2 /2m p h t- 

The system is driven out of equilibrium by its coupling to separate pumping and decay 
baths, so that the full Hamiltonian is given by H = H sys + H^^ p + H^?^. The contribution 

of the pumping bath is £™ = £ r ^A, n + %D j>n + H.c.) + £ >']j n).,,!),,, - & jtf fi jt J. 

j,n j,n 

The fermionic operators D^ n ,Cj n describe the pumping bath modes, and Tj jTl is the coupling 
strength. Similarly, the contribution of the decay bath is H^^ y = £ (p,p z (^p^p.Pz + H.c.) + 

P,Pz 

^2 ^p,?,— pj> z — pj>«> with H PiPz describing bulk photon modes. Each confined photon mode 

p couples to a separate set of bulk photon modes with various values p z , corresponding to 
conservation of in-plane momentum in the coupling between cavity and bulk photon modes. 



Path-Integral Formulation 

Following Ref. |14j . we construct the non-equilibrium generating functional Z as a coher- 
ent state path integral over field^] defined on the closed-time-path contour, C. For concise- 
ness, we arrange the fermionic fields into a Nambu vector A = (d, c) T . Formally, the par- 
tition function is thus: Z = J JT p pz ■ n T>[ty p , Aj, Cj, n , Dj^ n , E PtPz ]e lS , where the actiorj^] S = 
J c dt{^p{t)Kj(t)Cj, n {t)D^ n {t)^ pz ^ in terms of co- 

herent states specified by the fields *Sf p , Aj, Cj >n , Dj, m ^p,p z on * ne Keldysh time contour C. In 
what follows, we consider fields in the classical and quantum rather than forward and backward 
basis. 

1 In keeping with the convention of Ref. [13], we also refer to field amplitudes defined at discrete momenta, 
such as *P P , as fields. We note that Z is necessarily a functional integral, as we must account for a continuum of 
paths taken by $ ', A, C, D, S (and their complex conjugates) as functions of the continuous time variable t. 

2 When evaluating things we tend to take the continuum limit over p, making the partial time-derivative more 
convenient and appropriate [14] , 
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Treatment of Environment 

For the bath Hamiltonian given above, the action S contains only terms linear or quadratic 
in the bath fields Cj, n , Dj )n ,E, PiPz and their conjugates. Thus, the integral over these fields is 
Gaussian, and can straightforwardly be evaluated analytically. For the decay bath one thus 
finds: 

'oftfpCf). (2) 
- — p Pz 

where is a Pauli matrix in the (Keldysh) space of quantum and classical fields. By defi- 
nition, the bath has a large number of modes, and these modes thermalise rapidly compared 
to system-bath interactions. Hence we may take the bath occupation functions as fixed, and 
then allow the system distribution to be self-consistently determined. The Green's function 

of a free bosonic field is [(ihdt — hujpjaf]- 1 = (^t^_£/) Dp g where we have written 

P = (PiPz)- In frequency space the retarded and Keldysh Green's functions are given by 
D*(u) = [hco - ?u4 + id}' 1 , Dg(uj) = (-2m)F s (u;)S(u; - w£). Here F S (E) = 1 + 2n E (E) where 
n E (E) is the occupation function for the bath modes, which can have any form. For our situation, 
the decay bath is empty. 

Before proceeding further, we make some simplifying assumptions about the baths. We 
assume the bath frequencies Wp iPz form a dense spectrum, and the coupling constant ( PlPz = 
C( w p,p z ) is a smooth function. We may then replace summation over bath modes by integration. 
Then, taking the bath density of states and to be frequency independent, we have: 

O = £ *-Enc«> ( I 2t JS hw) ) »pM- <»> 

We follow an analogous procedure for the pumping baths (see [19] for details). 
Integration over Fermionic Fields 

After integrating over the bath degrees of freedom the full action S is: 



odecay 
'-'bath 



dtd* / £*p(*KE£ p . (ihd t 



//_: 



dtdt' 



3 P 

ihdt — 



where D 



(0),P 



(t, t') 



ihdt — hujp + in c 2in c F-={t — t') 



(4) 



To specify the exciton Green's function Gj we first introduce the abbreviations X c \(t) = gj^ Vfi i{t) / V% 
and Aq(t) = EpPj^W/v^ so that: 



a J 1 



( -Aq(t) i%d t -€j-ij x -A c i(t) ^ 

-A*(i) -A*j(i) ihdt + tlx 

ihdt-ej+i^x -Ad(t) 2ij x F D (ihd t ) -A q (t) 

\ -A*(t) ihdt + ej + vyx -A*(t) 2i lx F c (ihd t ) ) 



(5) 



where Fc\d(E) = 1 — 2nc t D(E) with tlc,d(E) the pumping bath occupation functions. As the 
occupation functions of all baths appear in this action, they compete to set the occupation 
function of the polaritons. This non-equilibrium action thus combines strong exciton-photon 
coupling with the effects of dissipation due to the open nature of the system. The action is 
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quadratic also in the fermionic fields Aj, so we can also integrate over these fields to get the 
effective action for the photon field alone: 

S = -i^Tr {inGj 1 } + / Hdidt' £ %(t)D7* (t, t')* p (t>). (6) 

j JJ-oo p 

As yet, we have made no assumption about what form \£ p (i) takes, however, since TrjlnGT 1 } 
involves ^ p (t), this effective action is nonlinear, so to proceed further analytically, some expan- 
sion or approximation scheme is required. Section [3] therefore discusses the mean- field theory of 
this model, and how it relates to laser theory as well as equilibrium results. 



3 Mean-Field Condition for a Coherent State 

The mean-field theory of the non-equilibrium system describes a self-consistent steady state, 
which may be found by evaluating the saddle point of S with respect to photon field, 5S/5^/ p d = 
SS/S^f* = 0. The first equation is satisfied if the quantum component vanishes, ^ p<g = 0. For 
the classical component, we write VPp cl = V^^p, so cfr p corresponds to the expectation of photon 
annihilation. If condensed, we consider the ansatz <fr p = (fro exp(— ifJ,st)8 Pt o, controlled by the 
parameters (fro, [is- For this ansatz to satisfy the saddle point equation, one requires 

(HuJo -us- *«c)^0 = 2 / 2^ G c)d^- ( 7 ) 



Putting the G^ d component of Eq. (|5| into Eq. (j7j) and defining E 2 - = (ej — ^ s /2) 2 + gjcfr^ we 
have the saddle point (mean-field) equation: 



(fko ~ fJ-s ~ iK c )(fro = ^29j(frolx 

j 

x f dv [F D (u) + F c {y)]v + [Fpjy) - FciyW, - fi s /2 + i lx ) 
J ^ [{ V -E j f + 1 l\[{v + E,Y + 1 l] ■ 1 > 

As noted above, the pumping bath occupations are imposed by choice, and we choose these 
to model a thermalised reservoir of high energy excitons, with a population set by the strength 
of pumping. In order to obey on average the constraint that we consider two-level systems, 
we take nc(v) + noty) = 1- Introducing parameters hb,Pb to describe the occupation and 
temperature of this exciton reservoir, we thus define: 



Fc,d{v) = tanh 



Pb ( , VB - Ms 
2 



(9) 



([is appears here via a gauge transform required to remove explicit time dependence from the 
effective action) . If there were no exciton-photon coupling the excitonic two- level systems would 
be thermally occupied, i.e. {d^d- — CjCj) = — tanh[/?s(ej — /is/2 + /is/2)/2]. 

As anticipated above, Eq. ([§]) is rather general, encompassing limits that correspond both to 
the equilibrium gap-equation for our model polariton system (discussed in Section [3J , as well as 
being capable of recovering the standard laser limit (discussed in Section [3J. In addition, if one 
extends this approach to slowly varying condensates, then as discussed in Section [3j one may 
make contact with the complex Gross-Pitaevskii approach. 
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Equilibrium Limit of Mean-Field Theory 

The simplest limit of the self-consistency equation, Eq. ([8]), is the thermal equilibrium limit, 
which corresponds to taking 7 x ,k c — > 0. In taking this limit, it is necessary to send k c — > 
first and then 7 X — > 0. This is because the self-consistency equation contained only the 
coupling of coherent photons to the decay bath, hence the decay bath cannot impose a non- 
trivial distribution on the system, while the pumping bath can. In order to satisfy Eq. ^ 
with k c = 0, the imaginary part of the right hand side must vanish. The most general way to 
achieve this is to set Fd{v) = Fc{f), which, considering Eq. Q, implies n$ = ^b- That is, in 
the absence of decay, one has chemical equilibrium between the pumping bath and the system. 
With fj>s = HB, the remaining part of Eq. <^ becomes: 

2i f dv 2 tanh (/3bi//2) v 

( "^ = V ^ ^W^WT^TWT^y (10) 

In the limit of small ~f x , one may use that lim 7a ._>.o ^7x/[( u ~ Ej) 2 + lx\ = 27r<5(V — Ej) to find 
(hiOQ — [ib)4>o = Yljidj&o/tEj) tanh (j3BEj/2) . This is the equilibrium mean-field theory [ToT - HT] 
of the system Hamiltonian introduced above. 



High Temperature Limit of Mean-Field Theory — Laser 

An alternative limit to thermal equilibrium is that of a simple laser. This limit too can be 
recovered from Eq. Q, in this case by taking Fcd^) to be frequency independent. This 
frequency independence can be recovered from Eq. ^ in the limit T — > oo, while keeping 
fiB oc T in order that the bath population remains fixed. Another interpretation of this is 
that infinite temperature corresponds to white noise, i.e. a Markovian approximation, where the 
occupation of the bath modes is frequency independent. [In contrast, Eq. ([8]) has a flat density 
of states of the bath, but a non-Markovian, i.e. frequency dependent, occupation]. 

Taking Fc^d{ 1 ') to be frequency independent, the integral in Eq. <JsJ) can then be simply 
evaluated by contour integration to give: 

(hu -pa- z. c )0o = Y,9 2 MFd ~ Fc f j ~^Jl^ X . (11) 

The term on the right hand side, describing the two-level system polarisation, is proportional to 
the bath inversion N$ = (tid — nc) = —{Fd — Fq)/2. In the limit cpo — > 0, this equation recovers 
the standard threshold condition for a laser [18j . This is clear if one restricts to gj = g,€j = e 
so the sum is replaced by a factor n, and one assumes resonance, 2e = Tiujq = /i, which yields: 
2K c j x /g 2 = uNq =total inversion. 



Low Density Limit: Complex Gross— Pitaevskii Equation 

Equation ^ is written for a uniform steady state, but in many cases, it is interesting to allow 
for solutions that vary slowly in time and space. To do this rigorously requires some care, 
but the basic idea can be described simply: One may start by writing Eq. ^ in the form 
(fj-s + iKc ~ fous)4> = x[0)/ i s , ]'^> where xl&i^s] is a nonlinear complex susceptibility. If one 
then separates the fast and slow time dependence ^(r, t) = 0(r, t)e~ tulot , one may write (ifhdt + 
ik c — [V(r) — h 2 V 2 /2m])(j)(r, t) = x[ < / ) ( r ) ^)]0( r > t), having introduced an external potential V(r). 
Then, by making a gradient and Taylor expansion of the nonlinear complex susceptibility xW\i 
one is naturally led to a complex Gross-Pitaevskii equation: ifidttp = (—h 2 V 2 /2m* + V c g(r) + 
U\(p\ 2 + i[7net(^B) — Tl'Al 2 ])'/') where T represents a nonlinearity of the imaginary part of the 
susceptibility. The dynamics of the excitons are responsible for producing an effective polariton 
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mass and effective potential. In some cases, it may also be important to consider the dynamics 
of the reservoir excitons more carefully, by introducing an extra degree of freedom to describe 
them. 



4 Applications: Fluctuations and Instability Towards BEC 

So far, we have discussed only the mean-field properties of the non-equilibrium polariton con- 
densate. We next consider fluctuations about this mean field, and in particular the photon 
Green's function. This is important for several reasons. Firstly, knowledge of the fluctuations 
determines whether a state is stable (i.e. do fluctuations grow or decay in time). Secondly, the 
photon Green's function describes the fluctuation contribution to physical observables such as 
luminescence and absorption spectra. 

To determine the photon Green's function, one may start from Eq. ([6]), and expand ^ = 
to second order in ip. The inverse photon Green's function has two parts, one from the bare 
photon action [given in Eq. Q], and one from expanding the trace over excitons. To determine 
the exciton part one may write Gj 1 = (G^ p ) -1 + 8GJ 1 , where Gj 3 is the fermionic Green's 
function including the mean- field photon field (f>, and 5GJ 1 is the photon fluctuation part given 
by: SGj 1 = —gj(ipqCr^ + ip q a^)a^/V2 — gj^^a^f + il) c \cr+ /v 7 ^ where a dc are Pauli matrices 
in the space of fermionic fields c, d. The action then depends on 



Tr {inGj 1 } = Tr jln [(Gf + GfSGj 1 - ^GfSGfGfSGfl 



(12) 



The last term gives a contribution quadratic in ip, which contributes to the inverse photon 
Green's function. 

When considering the condensed state, it is necessary to allow for anomalous correlations. 
This requires writing the Green's function in a (Nambu) vector space of (ipki^—k)' combined 
with the ± space due to the Keldysh/retarded/ advanced structure. Thus, in the condensed case, 
the Green's function is a 4 x 4 matrix, while when non-condensed it is only a 2 x 2 matrix. We 
begin by considering fluctuations in the normal state, and the nature of the instability to the 
condensate, and then in Section [4] briefly discuss fluctuations in the condensed state. 



Normal State Green's Functions and BEC Instability 

In the normal state, the spectrum and its occupation are determined by three real functions, the 
real and imaginary parts of the inverse retarded Green's function [D R (p, w)] _1 = A(p, u)+iB(uj) 
and the inverse Keldysh Green's function [D~ l (ijj)\ K = iC(uj). These functions can be read off 
from the fluctuation action. Once these functions are known, one may invert these expressions 
to find the retarded and Keldysh Green's functions, and thus determine the density of states 
p(p, a;) = — 2 Im[D R (p, lo)] and occupation of the modes 2n^[u) + 1 = iD K (uj)/ p(ui): 



2B(oj) 



CM 

2B(u) 



1 



(13) 



In terms of these, physical observables can be found, such as the luminescence C(p,uj) = 
p{p,uj)n^{oj). The roles of A(p,u), B(lj), C{oj) can be understood by considering the contribu- 
tion from the bare photon action. In this case A(p,uj) = Tloj — fiijo p determines the locations of 
the normal modes, B(uj) = k c gives the linewidth of these modes and C(uj) = 2K c [2ra=(u;) + 1] 
describes their occupation. Including the effect of the excitons, the zeros of A(p, u) now describe 
polaritons, rather than bare photons, and in addition B{u) is no longer constant, hence it plays 
a second role: If B(u) vanishes at some Twj = p e ff then this causes the occupation rif(fj, e e) to 
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Figure 1: Trajectories of zeros of A(uj, p = 0),i?(w) (i.e. normal modes £o and 
effective chemical potential fj, e g) for: (a) equilibrium weakly interacting dilute Bose 
gas (WIDBG); (b) equilibrium polariton condensate; (c) non-equilibrium polariton 
condensate; (d) Maxwell-Bloch laser. 




diverge. However, as long as A(p, fj, e g/h) does not vanish, the density of states will be zero at 
fj> e ff, so the luminescence will remain finite. 

A diverging occupation and vanishing density of states is exactly what would, in equilibrium, 
occur at the chemical potential, hence the identification ^ e g. Although the non-equilibrium 
system may be far from thermal, the emergence of zeros of B(oj) thus still describes an effective 
chemical potential. To ensure V p ^4(p, /J-eg/h) ^ 0, it is necessary that /U e ff is below any of the 
zeros of A(p,oj), i.e. the chemical potential is below all the polariton modes, as expected in the 
normal state. We next discuss the instability as fj, e g approaches the bottom of the polariton 
spectrum. If /x e g- is near the bottom of the polariton spectrum, we can expand A(p,u), B(oj) 
near their simultaneous zero, i.e. A(p,u) = a(hoj — S tP ), and B(u) = f3(hoj — /i e ff)- One may then 

(D R ) (D R ) 

find where the actual complex poles wj, of the retarded Green's function occur: fiuj^ = 
[(° 2 Cp + /3 2 /^ e fr) + ^o/3(/U fr — £p)][ a2 + P 2 ] ■ These poles determine the response to a small 
perturbation; thus, for perturbations to decay the poles must have a negative imaginary part. 
Hence, if ^ e fT > £p) then perturbations at that p grow and the normal state is unstable. One 
may also show that the point where /x e g- = £o> coincides with the first point where it is possible 
to satisfy the mean-field equation, Eq. ^ with </> = and fj-s = /^eff = Co- 
One can now understand the behaviour of the non-equilibrium system as pumping (and hence 
/ie) increases: At very weak pumping, where \ib is large and negative, B(oj) is always positive 
(i.e. decay dominates over gain), and no fj, e ff exists. As [ib increases, a region of negative B(uj) 
develops, and the boundaries of this region define fi e g as discussed above. As long as fi e g < £o the 
normal state remains stable. At the critical pumping power, /i e ff then reaches the lower polariton 
mode at p = 0, the normal state becomes marginally stable, and the mean-field equation can 
be satisfied. Beyond this point, the normal state would be unstable, but the condensed solution 
is now possible (and can be shown to be stable). 

Fig. [l|c) shows the evolution of /x e fr and £ for the non-equilibrium polariton condensate. For 
comparison, Fig. [l](b) shows the behaviour of the system Hamiltonian in thermal equilibrium 
and Fig. []Ja) that of weakly interacting and dilute Bose gas. One may note that despite the 
absence of a thermal distribution in Fig. ^c), the scenarios of normal state instability in these 
two figures are very similar, and would remain the same as long as the distribution function 
develops a divergence while the polariton system remains in strong coupling. In the next section, 
we discuss a case where the instability is somewhat different, that of the simple laser discussed 
in Section [3] and shown in Fig. [l](d) . 

Normal-State Instability for a Simple Laser 

As in Section [3j one may contrast the behaviour of the non-equilibrium condensate to that of a 
simple laser described by the Maxwell-Bloch equations [18] . which corresponds to the high tern- 
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perature (white noise) limit of our problem. As the above analysis concerns the retarded Green's 
function, we must define this function for the Maxwell-Bloch equations. The retarded Green's 
function describes the linear response of the system to an applied field, and so if one introduces 
a field Fe~ lult coupled to the photon field, then one has by definition ip(t) = iD R (u)Fe~ lult . One 
then finds: 

[Z^Hr 1 = A(lj) + iB{u) = hoj- hco + in c + V 9 j N ° - . (14) 

3 

For this form, B(uS) can only become negative if Nq is sufficiently large. Restricting as in Section 
[3] to gj = g, ej = e = Tiujq/2 one finds the requirement for gain is g 2 nNo > 2k c <j x , which is again 
the laser threshold condition. In this same restricted case, the zeros of A{uj) behave as follows: A 
solution £ = always exists, and if Nq < —^ 2 /g 2 n, an extra pair of roots exist. The evolution 
of these zeros is shown in Fig. [ijd). One may note that for the Maxwell-Bloch equations, strong- 
coupling (i.e. splitting of the modes £) collapses before condensation (lasing) occurs, while for 
the polariton condensate, condensation occurs while still strongly coupled. 



Fluctuations of the Condensed System 

As noted earlier, if condensed, the Green's function is a 4 x 4 matrix, so the derivation of the 
spectrum becomes more complicated, however the essential features can be explained by general 
arguments. The following discussion is thus based on symmetry arguments. (The full derivation 
matches these results [191). The form of the inverse retarded Green's function is constrained by 
the following requirements: there must be symmetry under p — > — p; the modes must in general 
have a finite linewidth; however at p — > 0, there must be a mode with vanishing frequency 
and vanishing linewidth corresponding to global phase rotations of the condensate. These three 
considerations determine the leading order behaviour of D R (p,uj) for small a;, p. Using these 
ideas, one may then write: 

(P ' U) = detQ^p) = uj 2 + 2icjx - c 2 p 2 ' (15) 

The parameters x, c describe the linewidth and sound velocity. 

^From this form of D R (p, to), one finds the poles are given by wj/ 5 ^ = — ix±i\J x 2 — c 2 p 2 . At 
long wavelengths, these are diffusive (only an imaginary part exists), and only above a critical 
momentum does a real part emerge. Given the generality of the argument leading to this 
result, it is unsurprising to find the same structure emerges from other approaches, see e.g. [20] . 
Similar results also occur for the case of a parametrically pumped polariton system |21| . The 
absence of a linear dispersion of energy vs momentum in the condensed state affects some aspects 
of superfluidity in this non-equilibrium system, however there are also aspects of superfluid 
behaviour that survive | 



5 Connection to Other Approaches 

The language of the Keldysh path integral, and the Keldysh Green's functions provide a natural 
bridge to many other approaches that have been used to model non-equilibrium polariton con- 
densation. We have already discussed above the connection between the mean-field theory, i.e. 
the saddle-point of the Keldysh action, and the complex Gross-Pitaevskii equation [20l l24"t [25] . 
In order to go beyond mean-field theory, the approach discussed in this chapter makes use of 
Keldysh/retarded/advanced Green's functions to describe both the occupation of a mode, and 
the density of states. These Green's functions can naturally be related to the one particle density 
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matrix p(r, r', t) = (^(r,t)^(t,r')) = i(D K -D R +D A )(r, r',t, t)/2. Direct time evolution of the 
one particle density matrix has been used to treat polaritons in zero-[26j and one-dimensional 
geometries [27], and stochastic methods for simulating density matrix evolution have been used 
to describe a number of properties of polariton condensates, see e.g. [251 Another stochastic 
approach used for exciton condensation is the Heisenberg-Langevin equations |30j . Such an 
approach again connects naturally to the Keldysh formalism, with the Keldysh self energy due 
to the bath corresponding directly to the noise correlator of the Heisenberg-Langevin approach, 
and the retarded self energy corresponding to the dissipation term. Finally, there is also a clear 
connection between the Keldysh Green's functions and the quantum Boltzmann equation (see 
e -g- jSH [32] for a derivation of the quantum Boltzmann equation from the equations of motion 
for the Keldysh Green's functions). There have been many works using the Boltzmann equation 
to model kinetics of polariton condensation [33H40] . By considering how the quantum Boltz- 
mann equation arises from the Keldysh Green's function, one may note that in order to correctly 
describe the coherence properties of the condensed state, one must include anomalous retarded 
self energies, modifying the polariton spectrum. 
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